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Abstract 

These lectures give an introduction to the methods of conformal 
field theory as applied to deriving certain results in two-dimensional 
critical percolation: namely the probability that there exists at least 
one cluster connecting two disjoint segments of the boundary of a 
simply connected region; and the mean number of such clusters. No 
previous familiarity with conformal field theory is assumed, but in the 
course of the argument many of its important concepts are introduced 
in as simple a manner as possible. A brief account is also given of 
some recent alternative approaches to deriving these kinds of result. 



*Lectures delivered in "New Trends of Mathematical Physics and Probability 
Theory", Chuo University, Bunkyo-ku, Tokyo, March 5-6, 2001. 



1 Introduction. 



The percolation problem has for many years been of great interest to theo- 
retical physicists and mathematicians, in part because it is so simply stated 
yet so full of fascinating results. It embodies many of the important concepts 
of critical phenomena, yet is purely geometrical in nature. 

Percolation studies the clustering properties of identical objects which are 
randomly and uniformly distributed through space. In lattice bond percola- 
tion, the links of a regular lattice, of edge length a, are either open or closed. 
In the simplest version of the model, the open bonds are independently dis- 
tributed with a probability p for each to be open (and 1 — p to be closed.) 
In site percolation, the bonds are all assumed to be open, but now each site 
is open with probability p. In both cases we study the statistical properties 
of clusters of neighbouring open bonds and sites. When p is small, the mean 
cluster size is also small, but, in more than one dimension, there is a critical 
value pc of p, called the percolation threshold, at which the mean cluster size 
diverges. For p > Pc there is a finite probability that a given site belongs to 
an infinitely large cluster.]^ 

In these lectures, we shall be concerned with properties of the contin- 
uum limit of percolation. This may be defined as follows (we consider two 
dimensions from now on): consider a finite region TZ of the plane, bounded 
by a curve F. Consider a percolation problem on a sequence of lattices Ca 
covering TZ, constructed in such a way that the lattice spacing a keeping 
the size of TZ fixed. Obviously when we do this, such quantities as the total 
number of clusters in TZ will diverge as a — 0, so we need to identify some 
suitable quantities which might have a finite limit. An example is afforded by 
the crossing probabilities. Suppose for simplicity that 7Z is simply connected, 

^ There are also so-called continuum versions of percolation, for example the clusters 
formed by hard spheres of a given radius a which are distributed independently so that 
they may overlap. In all cases, however, a finite microscopic length a is necessary to define 
the notion of clustering. 
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Figure 1: A crossing cluster from 71 (C1C2) to 72 (C3C4). 

and let 71 and 72 be two disjoint segments of F (see Fig. Then a crossing 
event is a configuration of bonds (or sites) on the lattice Ca such that there 
exists at least one cluster, wholly contained within TZ, containing both at 
least one point of 71 and of 72. Let P(7i,72;£a) be the probability of this 
event. Then the following statements are all conjectured to be true: 

• P(7i,72;£a) has a finite limit, denoted by -P(7i,72), as a ^ 0, which 
(under a broad class of conditions) is independent of the particular 
form of £, of the precise way in which C intersects the boundary F, 
and of whether the microscopic model is formulated as bond, site, or 
any other type of percolation, as long as there are only short-range 
correlations in the probability measure.0 

• P(7i,72) is invariant under transformations of TZ which are conformal 
in its interior (but not necessarily on its boundary F). 

• The Riemann mapping theorem allows us to conformally map the in- 
terior of TZ onto the interior of the unit disc \z\ < 1 of the complex 

■^Of course, this limit is interesting only at p = pc- For p < pc, linia^o ^ = 0, because 
all clusters are finite (in units of a), while for p > pc the infinite cluster always spans, so 
the limit is 1. 
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plane. Suppose that the ends of the segments are thereby mapped into 
the points (zi, Z2, -23, Z4) (see Fig. ^for the labeUing). Then P(7i, 72) is 
a function only of their cross-ratio 

"I - {Z, - Z,){Z2 - z,) ' 



and has the explicit form 



r(|) 
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^=^iy^^^2Fi(|,|;|;ry) (2) 
where 2F1 is the hypergeometric function. 

Moreover, if the random variable A'"c(7i,72; Co) denotes the total num- 
ber of distinct clusters which cross from 71 to 72, then the whole prob- 
ability distribution of Nc also has a finite limit as a — and is confor- 
mally invariant, depending only on rj. In particular the mean number 
of such crossing clusters is 
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(3) 



Equations are just two among many similar results which may been 
obtained using methods of conformal field theory developed by theoretical 
physicists fl], ||, H. At first sight, the appearance of formulas like these may 
seem quite mysterious to those used to thinking about percolation as a lattice 
problem. It is certainly true that the methods originally used to derive them 
are not, so far, mathematically rigorous.0 But there is no doubt that they 
are correct - (0) has been numerically tested to great precision in a num- 
ber of casesQ, and, moreover, once the existence of a conformally invariant 
continuum limit is accepted, the formulas follow from very classical mathe- 
matical methods. It is the purpose of these lectures to give some idea to a 
^S. Smirnovpl has recently given a proof of (p|) using other methods. 
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non-specialist audience of how these kinds of resuhs arise. Inevitably I shall 
not be able to cover all the details, but hopefully the lectures will provide a 
basis from which to explore the literature further [Q.^ For reasons of time I 
will focus on the derivation of the above formulae, which is mostly my own 
work, as well as mention some other recent alternative derivations. This is 
not to overlook the work of others on other important results in percolation, 
particularly that of Nienhuis, Duplantier, Saleur and others on the Coulomb 
gas approach, which is better suited to systems without boundaries. 

2 Percolation, the random cluster model, and 
the Potts model. 

The emphasis in this course will be on the analogy between percolation and 
conventional critical behaviour in spin systems. This is through a well-known 
mapping first discovered by Fortuin and Kastelyn[0. The Potts model is a 
generalisation of the Ising model in which the spins s(r) at each site of a 
lattice take the values {1,2, ... ,Q), where, initially, Q is an integer larger 
than 1. The energy is the sum over all nearest neighbour pairs (r',r") of 
sites (i.e. a sum over all bonds) of —JSs(r'),s{r")- Thus the partition function 
is 

Z = Tl exp(/3J ^ 5s{r'),s{r")) (4) 

Apart from an overall unimportant constant this may be rewritten as 

Z = Tr Yl ((1 - P) + pSs{r'),s{r")) (5) 

where p = 1 — e^^'^ . 

^The main text includes only those concepts and arguments needed to arrive at the final 
conclusion. Footnotes will indicate how these ideas fit within the more general framework 
of conformal field theory. 
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Now imagine expanding out the product. If there are B bonds on the 
lattice there will be 2^ terms in this expansion. Each term may be associated 
with a configuration in which each bond of the lattice is open (if we choose 
the term (x p), or closed (if we choose the term oc (1 — p))- Sites connected 
by open bonds form clusters, and the Kronecker deltas force the all the spins 
in a given cluster to be in the same state. When we trace over the spins, 
each cluster will have only one free spin, so will give a factor Q. Thus we 
can write Z as a sum over configurations C of open bonds: 

Z = ^p|c|(l-p)^HC|Q^('^) (6) 
c 

where N{C) is the number of distinct clusters in C. Note that in this form, 
at least for a finite lattice, Z is a polynomial in Q and therefore its definition 
may be extended to non-integer values of Q. 

The weights in (^ define the random cluster model. Of course, in perco- 
lation, each cluster is weighted with a factor 1, so it corresponds to Q = 1. In 
that case, the sum is simply over all possible configurations weighted by their 
probabilities, so Z{Q = 1) = 1. However, there is nontrivial information in 
the correlation functions. For example, the probability that sites ri and r2 
are in the same cluster is given by the limit as Q — ^ 1 of 

where a and b are any two different Potts states.^ From (^) we can also 
calculate quantities like the mean total number of clusters by differentiating 
with respect to Q: 

E[N] = {d/dQ)\Q=^ZiQ) (8) 

It is important establish how quantities like the crossing probabilities 

may be related to those in the Potts model. For this we need to bring in the 

^ This follows from a similar argument to that which gives the crossing probabilities in 
terms of partition functions, described later. 



5 



notion of boundary conditions. Consider a Potts model defined on a lattice 
L which covers the region 7^ of the plane, as described in the Introduction. 
We consider the boundary of L as being the set of sites which lie just outside 
7^ but are adjacent to sites within 7^. These boundary sites form a discrete 
approximation to the boundary F. On these boundary sites, there are two 
simple and natural boundary conditions we might impose on the Potts spins: 
/ree, which means that we sum freely over them in the partition function; 
and fixed^ in which case they are all fixed into some particular Potts state, 
say a. 

Notice that if we take the same boundary condition on the whole bound- 
ary, then in the limit Q — > 1 partition functions with either free or fixed 
boundary conditions become the same, Z = 1, since all spins must be in the 
same state! But we can get something nontrivial if we allow the boundary 
conditions to be different on different parts of the boundary. In particular, 
consider the geometry of Fig. |I], and suppose the boundary spins are fixed 
on the segments 71 and 72, and free on the rest of the boundary. By the 
permutation symmetry of the Potts states, there are two possible different 
cases: when the fixed states on 71 and 72 are the same, for example a; or 
when they are different, say a and b. Let us denote the partition functions 
in the two cases by Zaa{Q) and Zab{Q) respectively. Now each configuration 
C either has, or does not have, a cluster which spans between 71 and 72 (see 
Fig. Those which do not have a spanning cluster contribute to both par- 
tition functions (with the appropriate weights), but configurations which do 
cannot contribute to Zab{Q)i since the existence of the cluster would force 
spins on the two segments to be in the same state, which is not true by 
assumption. Thus we have a simple relation for the crossing probability, on 
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Figure 2: Crossing events contribute to Zaa but not to Zghi while non-crossing 
events contribute equally to both. 
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the lattice:^ 

P(7i, 72; C) = Urn {Zaa{Q) - Zab{Q)) (9) 

Remember that each of these is a polynomial in Q, so that Zab makes perfect 
sense at Q = 1, even though it might seem that we need at least two states 
to define it.[] 

With only slightly more effort we can relate the mean number of distinct 
crossing clusters E[Nc\ to partition functions. Consider once again the geom- 
etry of Fig. |1], but now subdivide the clusters into those which touch neither 
7i nor 72, those which touch 71 but not 72, those which touch 72 but not 
71, and finally the crossing clusters, which touch both 71 and 72. Denote the 
number of each such cluster in each configuration by A'^o, Nl, and Nc 
respectively (see Fig. |^.) Since clusters which touch a boundary where the 
Potts spins are fixed are counted with weight 1, while those which are free 
are counted with weight Q, it follows that 

Zf^ = {Q^^+Nl+nm^ Z,„ = (Q^») (10) 
Zaf = (Q^«+^°) Z;,= (Q^-+^«) (11) 

where / now denotes that the boundary spins are free on that portion of 
the boundary (they are always free on the complement of 71 U 72.) Then 
straightforward algebra shows that 

(iV,) = {d/dQ)\Q=,iZff + Z,, - Zf, - Z,f) = {d/dQ)\Q=,{ZffZ,jZf,Z,f) 

(12) 

where the last equality holds because all the partition functions equal 1 at 
g = 1. 

^Note that the crossing probabiUty in the random cluster model with Q 7^ 1 is not given 
by the (normalised) difference {Zaa — Zab)/Zaa, since this counts clusters which touch 71 
and/or 72 with weight 1, rather than Q. 

^We can think of the crossing probability as P{Nc > 1), where Nc is the number of 
crossing clusters. It is also possible to design more complicated boundary conditions which 
give P{Nc > n) for any integer n> 1. See 
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Figure 3: Different types of cluster, which touch both, one or none of 71 and 

72- 

It will be the aim in the rest of these lectures to derive explicit expressions 
for these partition functions, in the continuum limit, and for Q close to 1. 

3 The continuum limit of critical lattice mod- 
els. 

The Potts model in two dimensions is known to have a critical point, with 
a divergent correlation length, for < Q < 4. We now move onto less solid 
ground. What I am now going to assert is based on evidence from the exactly 
solved Ising model Q = 2, as well as the analysis of renormalised perturbation 
theory near the upper critical dimension. ||] 

At the critical point, these theories are believed to be scale invariant. 
What this means for correlation functions is the following. Consider for ex- 
ample, the correlations (s(ri)s(r2) . . . s(r„)) of the local lattice magnetisation 
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in the Ising model. Since s{r) = ±1, this quantity has no dimensions. At the 
critical point, in the limit where the lattice spacing a — > with the points rj 
kept fixed, it behaves like a"-^ times a function of ri, . . . , r„, where x is a pure 
number (equal to | for the Ising model magnetisation.) This means that we 
can define a scaling operator (f){r) = a~^s{r) such that the limit a — > of the 
correlation functions (0(ri)0(r2) . . . (p^rn)) exists. The pure number x, which 
we should denote called the scaling dimension of (f). 

Moreover, as long as the original lattice model has sufficient symmetry un- 
der finite rotations, the limit is invariant under infinitesimal rotations. Thus, 
far away from any boundaries, the two-point correlation function (0(ri)0(r2)) 
depends on the separation |ri — only, and, on dimensional grounds, must 
therefore have the formQ 



The above statement is strictly true only when the points rj are distinct, 
in the limit a ^ 0. That is, they are separated by an infinite number of 
lattice spacings. When their relative distances are kept fixed in units of a, 
other scaling operators may arise. For example, the product s{r[)s{r'() on 
neighbouring sites may be thought of as being proportional to the local energy 
density in the Ising model. Correlation functions of this quantity behave in a 
similar way to those of the magnetisation, but they define a different scaling 
operator, with a different scaling dimension. In general there is an infinite 
number of scaling operators 0fc(r), each with their own scaling dimensions Xk- 
Arbitrary local products 5*^ of spins which are separated by distances 0(a), 
are given asymptotically as linear combinations of these scaling operators: 



k 

where the dimensionless coefficients Ajk are of order unity, and correlations 
of the scaling operators all have a limit as a — > 0. These are rotationally 

^In the case of the magnetisation, 2x^ is conventionally denoted by {d~ 2 + i]). 



{(f){ri)(j){r2)) = const. |ri — r2 



2x0 



(13) 




(14) 
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and translationally invariant (far from a boundary), and also scale covariant 
under a scale transformation r ^ r' = h~^r 

(0i(r;)02(r^) . ..<i>nK)) = 6^^"M0i(ri)02(r2) • • .0n(r„)) (15) 

In ([l^ the coefficients depend on the particular lattice and so on, but the 
scaling dimensions are universal. For example, in the Potts model they 
depend only on Q and the dimensionality of the system, assuming that the 
interactions are short-ranged.^ 

3.1 Boundary operators. 

The above statements about the existence of the continuum limit remain valid 
in the presence of boundaries. Of course the boundary now breaks global 
rotational and translational invariance. Taking for simplicity the boundary 
to run along the real axis, so that the system occupies the upper half-plane 
y > 0, then the two-point function of a scaling operator has the general 
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(0(a;i,l/i)0(x2,l/2)) = (1/1Z/2) """^F {\xi - X2\/yi,\xi - X2\/y2) (16) 

where F is a scaling function. This is valid for all yi and y2 strictly positive, 

in the limit a —>■ 0. F behaves in such a way as to recover the bulk result 

((xi — X2)^ + (2/i — 2/2)^)"^"^ as both yi, 7/2 — ^ 00, but it turns out that the limit 

^ One of the important properties of the scahng operators is that they are complete in 
the foUowing sense. If, in the above correlation function, we take the separation |ri — r2 1 to 
be much smaller than the distances |ri — r^j to the other points with fc > 3, this correlation 
function may be expressed as a sum over correlation functions with only a single operator 
at the mid-point, say, of ri andr2: (0i(ri)02(r-2) ■ • ■) ^J2kCi2k\ri-r2\'^''~'^^~'^^(j)k{{ri + 
r2)/2) . . .) where the dependence on ri — r2 is dictated by scale covariance and rotational 
invariance (we have assumed that all the operators are scalar under rotations, which is 
a slight over-simplification). The important thing is that the coefficients Cijk are also 
universal, and in particular independent of the other operators hidden in the .... We can 
thus remove the (...). The result is called the operator product expansion (OPE). 
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as yi and/or ?/2 — ^ is singular. In that limit we find instead the behaviour 

k 

where the Bk are (universal) constants and the set of boundary scaling 

dimensions. Thus there is a separate set of scaling operators (pk defined on 
the boundary y = 0, which have scaling dimensions which are in general 
different from those in the bulk.|^ It turns out that the set of such operators 
depends not only on the bulk universality class, but also on the nature of the 
boundary condition, for example whether the lattice spins are free or fixed. 

3.2 Boundary condition changing operators. 

Once the idea of boundary scaling operators was understood, it was realised JIT 
that, at least in two dimensions when the boundary is one-dimensional, there 
is another set of objects which should possess similar properties under scale 
transformations. These are points on the boundary where the boundary 
condition changes from one type to another. Let us illustrate this with a 
relevant example from the Potts model. Consider a Potts model in the half- 
disc r < R, y > 0, see Fig. ^j. Compare the case when (a) the boundary 
conditions on the Potts spins are free (/) at every point of the boundary 
except for the interval {zi, z^) of the real axis, where they are fixed to state 
a, to the case (b) when they are free everywhere on the boundary. Denote 
the corresponding partition functions by Zaj and Z^. Then the ratio Zaf/Zf 
has a finite limit as ^ oo, |^ which therefore depends only on \zi — Z2\. 
We may think of this ratio as defining the two-point correlation function of 
boundary condition changing operators (bcc operators) with a well-defined 
scaling dimension: 

(0/|a(^l)0a|/(^2)) = ]\ma-^^'^{Zaf/Zf) CC \z, - Zsl"'"^ (18) 

^"^(P^) is equivalent to the bulk-boundary expansion (j){x,y) = Bky^^'*^^''4'k{x). 
^^This is not obvious, at criticality, and needs to be proved. 
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Figure 4: Potts model in the half-plane. The two-point function of the bcc 
operator counts clusters which touch the real axis in the interval {zi^ Z2). 

Within the Potts model, we expect the scaling dimension of this bcc operator 
to be universal, depending only on Q. More complicated sets of boundary 
conditions may be considered. For example, the partition function Zab in 
(P) which we shall need to compute the crossing probability is related to the 
correlation function {(j)f\a4'a\f4'f\b4'b\f)- 

3.3 Finite-size scaling of the transfer matrix. 

Consider now a lattice model defined on an infinitely long strip of width L, 

parametrised by the coordinates (m, v) with < f < L. We could consider 

periodic boundary conditions in the f-direction, but for our purposes it will 

be more useful to think about particular boundary conditions, labelled say 

by a and (3 (which do not have to be the same) on either edge of the strip. A 

convenient way of discussing this geometry is through the transfer matrix]^ 

Tap{L), which is a finite matrix whose rows and columns are labelled respec- 

^^I shall attempt to be consistent and denote true operators by a hat ("). Note that the 
scaling 'operators' introduced earlier are not true operators, but simply local densities, 
which commute with each other. 
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tively by the states of the spins in neighbouring rows {u,u + a), and whose 
elements are the Boltzmann weights for the two rows. If we take the strip to 
have finite length W in the w-direction, with periodic boundary conditions 
in that direction, the partition function is given by 

Z = Trf^/« = 5^Aj;^/" (19) 

n 

where the A„ are a complete set of eigenvalues of T (which of course depend 
on a, j3 and L.) In this picture the lattice spins s{u,v) themselves become 
matrices s{v) (diagonal in this basis), acting on the same space as does T. 
If the eigenstate corresponding to the largest eigenvalue of T is denoted by 
|0), then a two-point correlation function may be written 

{s{u2,V2)s{u,,v,)) = (0|s(i;2)(T/Ao)("^-"^)/«s(t;i)|0) (20) 

and so, by inserting a complete set of eigenstates of T, we sec that it decays 
along the strip as a sum of terms of the form (An/Ao)*-"^^"^-*/" such that 
(ri|s|0) 7^ 0. Similar equations hold for the correlation function of any local 
product Sj{v) of lattice spins. 

How do these quantities behave in the continuum limit a ^ 0, with L 
fixed? The space of states becomes infinite-dimensional, but we shall not 
worry too much about its precise structure. The scaling 'operators' (j){u,v) 
become true local operators 0(t') acting on this space. It is useful to write 
the transfer matrix itself as 

fap{L) = eM-aHap{L)) (21) 

where H has eigenvalues En, so that correlation functions decay as a sum of 
terms e~^^"~^''^^"2-wi)_ ^i^q existence of the continuum limit implies that the 
gaps En — Eo have a limit as a ^ 0. Scaling then implies that, at the critical 
point, each of these gaps must be proportional to L^^. As we shall soon 
see, conformal invariance relates the constants of proportionality to scaling 
dimensions of operators. 
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4 From scale invariance to conformal invari- 
ance. 



Let us recall the equation (|l^) for the behaviour of a general correlation 
function under a scale transformation r ^ r' = b^^r: 

{Mr'i)Mr2) ■■■) = b^>=^HMri)Mr2) ■ ■ ■) (22) 

One way to understand this equation is from the existence of the continuum 
limit a ^ 0, and the absence of any other length scale at the critical point: 
if this limit does exist, it should not matter whether we rescale the lattice 
spacing in the underlying lattice model by a factor 6: a — ba, in the sense 
that the measure of those degrees of freedom which survive the continuum 
limit should be unchanged. But we can think of this rescaling (which is a 
kind of renormalisation group transformation) as keeping the lattice spacing 
fixed, and rescaling r ^ r' = b^^r. 

Now we make the following bold generalisation: since the measure is 
ultimately a product over local Boltzmann weights, we may equally well 
make a rescaling a — * b{r)a where the factor b now varies smoothly with r 
(on a scale ^ a.) So we generalise (E2p to 



{Mr'i)Mr2) ■ ■ ■) =l[KrkrHMri)Mr2) ■ ■ ■) (23) 

k 

where b{r) = \dr/dr'\, the jacobian of the transformation. 

What are the allowed transformations r — > r'? Locally, they must corre- 
spond to a mixture of a scale transformation together with a possible rotation 
and translation. Such transformations are called conformal^ An example 
is shown in Fig. So (^) describes the conformal covariance of an ar- 

Strictly speaking, conformal transformations operate on the metric rather than the 
coordinates, but this distinction is not important here. 

^"^It is not possible to have a non-trivial conformal transformation without some element 
of local rotation as well. This means that anisotropic critical points, which occur, for 
example, in directed percolation or at a Lifshitz point, are not conformally invariant. 
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Figure 5: Example of a conformal transformation, locally equivalent to a 
scale transformation plus a rotation. 



bitrary correlation function.[^ Once again, it is not rigorously founded, but 
rather abstracted from exactly solved cases like the Ising model, and from 
the structure of renormalised perturbation theory. It is generally believed 
that the continuum limit of any model which is invariant under scale trans- 
formations, rotations and translations, and has short-ranged interactions, is 
also conformally invariant. We are going to assume that it is valid for the 
continuum limit of the critical Potts model with < Q < 4. 

In two dimensions, conformal invariance is particularly powerful, because 

-"^^This is written for the case of operators which transform as scalars under rotations. A 
deeper analysis also shows that (^) cannot be true for all scaling operators: for example 
it is easy to show that if it is true for it cannot in general be true for V^0. How- 
ever, this analysis shows that it is true for a subset of operators called primary, and the 
transformation laws for all other operators may be derived from these. 
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if we label the points of the plane by a complex number z = x + iy, any 
analytic function z z' = f{z) defines a conformal transformation, at least 
at points where f'{z) ^ 0. 

4.1 Some simple consequences. 

Let us consider for definiteness correlations of boundary operators, since it is 
these we need for our problem. They may be ordinary limits of bulk operators 
at the boundary, or they may be bcc operators. Just as scale covariance 
fixes the form of the two-point functions, conformal invariance gives further 
information about the higher-point functions. If we consider the upper half 
plane, then the set of conformal transformations which preserve this are the 
real Mobius transformations 

^ = '^^±}_ ad-bc^ 0, (24) 

cz + a 

where a, b, c and d are real (this is to preserve the real axis.) Given a 4-point 
function {(f){zi)(f){z2)4>{z3)(j){z4)) , with the zj real, it is therefore possible in 
general to choose the parameters in (pif ) such that Z2, z^ and z^ are mapped 
into three pre-assigned points, say, (0, oo,l), corresponding to z' = {z — 
^2) (-2^4 — ^3)/(^ ~ -23) (-24 ~ -22)- However there is always an invariant of such 
a transformation: it is the cross-ratio (anharmonic ratio) 

rj = {zi2Z'ii/ Zi'iZ2i) (25) 

(where the notation Zki = Zk — zi has been introduced) so that, in this 
example, zi gets mapped into rj. We may now apply the transformation 
formula (p3D. The correlation function on the left is some (as yet) unknown 
function of r], and it is simply a matter of working out the jacobian. After 
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some algebra,^ 

(0(^00(^2)0(^3)0(^4)) = ( F{r^) (26) 

4.2 Application to crossing probabilities. 

We are already in a position where we can deduce part of the main result. We 
begin with the partition functions corresponding to the geometry of Fig. |l], 
in which the Potts spins on the segments 71 and 72 are fixed, and the rest 
are free. As we showed earlier, the crossing probability is given by the limit 
as Q ^ 1 of the partition functions Zaa and Zab- According to our analysis 
above these may be written, as a ^ 0, 

Zaa/Zf ~ a^-'(«)(0;|,(zi)0,|;(;22)0/|a(^3)0a|/(^4)) (27) 
Zab/Zf ~ a^^'(«)(0/|,(zi)0,|/(^2)0/|6(^3)06|/(^4)) (28) 

Here we have introduced the scaling dimension x{Q) of the (free | fixed) bcc 
operator, which will depend on Q, but not on a or b. Now the Riemann 
mapping theorem assures us that there exists a conformal mapping of the 
interior of the region TZ to that of the unit disc, and thence, by a simple 
Mobius transformation, to the upper half plane y > 0. Thus we may use (^) 
to relate the ratios of partition functions in the two geometries.P^ However, 
this will bring in non-trivial factors from the jacobian, evaluated at the points 
Zj, raised to the power x{Q), so that the ratio of partition functions is not, 
in general, conformally invariant. 

^^This is the simpler case when all four operators have the same scaling dimension. A 
more general result also holds. 

^^In doing this we must assume that the transformation is conformal also at the points 
Zj. This requires that the boundary T be difFerentiable at that point. There are interesting 
additional factors when the boundary has a corner at one of these points (as happens in the 
relevant case of the rectangle) . However all these factors are raised to a power proportional 
to x{Q), so are not important at Q = 1. See Ref. |l3| . 
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However, a remarkable thing happens at Q = 1. Consider the two-point 
function {(j)f\a{zi)(j)a\f{z2)) in the upper half plane geometry. As discussed 
earlier, this is proportional to the ratio of correlation functions Za/Zf, and, 
at the critical point, decays as \zi — 2:2 1"^"^*-'^-'. In the random cluster model, 
Za counts all clusters which touch the real axis in the interval {zi,Z2) with 
a weight 1, while counting all other clusters with weig ht Q (see Fig. |.) 
Denoting the total number of each type of cluster by A^i and N2 respectively, 
we have Za = {Q^^) and Zf = {Q^^'^^^). Trivially, then, both partition 
functions equal one at Q = 1, so that their ratio is independent of zi and 2:2- 
Thus 

a;(l) = (29) 

an innocent but remarkably interesting equation from the point of view of 
conformal field theory.0 From our perspective it means that the crossing 
probabilities in percolation (at least those which may be expressed in terms 
of 4>f\a) are conformally invariant, in the sense stated in the Introduction. 
Thus percolation has stronger properties under conformal transformations 
than do more general critical theories. This may be traced to the fact that 
its partition function is equal to 1. 

5 Scaling operators and states in conformal 
field theory. 

We have already discussed the continuum limit of the transfer matrix on 

strip of width L, and introduced the concept of a state space spanned by 

the eigenvectors of the operator H. Now we shall see that this has a very 

simple structure due to conformal invariance. Consider the upper half plane, 

with boundary conditions labelled by a along the real axis. Consider the 

^^Note that if we ask for the mean number (A^i) of clusters touching the interval (zi, Z2) 
we must first differentiate with respect to Q. The result is 2a:' (1) In |zi — Z2|. As we show 
later, the prefactor is V3/47r. 
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conformal mapping 



z = x + iy^w = u + iv = (L/tt) In (30) 

which maps the upper half plane into such a strip, with boundary conditions 
a on both edges. If we apply the transformation law (|23|) to the correlation 
function {(f){zi)(j){z2)) in the half-plane (where is some boundary operator 
and we assume that < zi < ^2), we find 



((^(Mi,0)<^(n2,0)) 



vr/L 



sinh((7r/L)|ui — M2I 



4> 

(31) 



For |mi — ^21 -C -Z^ this reduces to the same result as in the half-plane, but 
in the opposite limit it decays as exp(— (7rX(/,/L)|ui — M2I). This means that 
there must be an eigenstate of H with a gap En~Eo = ttx^/L. In fact it may 
be argued that to each boundary scaling operator there is an eigenstate. 
This may be identified by taking ui —00 (so that xi ^0), in which case 
only the state {n\ with the lowest value of En such that (?t,|0|O) 7^ will 
contribute. A less obvious statement is that to each eigenstate there is an 
operator. This can be seen by choosing the strip to be in some particular 
eigenstate \n) at u = ±00 and computing correlation functions on the strip 
with this state rather than the vacuum |0). On transforming these back to 
the half-plane, they turn out to be exactly as if we had an operator at the 
origin, with scaling dimension {L/7T){En — Eq). Thus we have one of the 
most important results of two-dimensional conformal field theory:]^ 

• There is a (1-1) correspondence between the (boundary) scaling oper- 
ators and the eigenstates of H in the strip. 



A more field-theoretic way to state this is that if H generates translations in u along the 
strip, then {L/-k){H — Eq) generates scale transformations in the half-plane. Eigenstates 
of the latter correspond to scaling operators. 
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In what follows, we shall normalise so that Eq{L) = O.Q This means 
that a boundary scaling operator (p, inserted near the origin, corresponds 
to an eigenstate with E = vrx^/L in the strip. We may extent this to bcc 
operators: if we now consider a strip with different boundary conditions a 
and (3 on either edge, this corresponds to an insertion of the bcc operator 
4>a\f3, and so the the lowest eigenvalue of Hai3 is nx^^^^/L. 

5.1 Descendent operators. 

Let us consider in more detail the structure of the set of scaling operators. It 
is useful first to think about infinitesimal conformal transformations (which 
preserve the upper half-plane). The simplest ones are the translations, z 
z + a_i (where a_i is an infinitesimal parameter), under which a scaling 
operator 0(0) transforms according to 

0(0) ^ 0(a_i) = 0(0) + a_i[L_i0](O) (32) 

where L_i0 is nothing but {d/dz)(j). There are also the scale transformations 
z ^ z + aoz, under which 

0(0) ^ (1 + ao)"^0(O) ~ 0(0) + ao[Lo0](O) (33) 

where Lo0 = x^0. The statement that the correlations of satisfy ( p3D is 
equivalent to assuming that for a general infinitesimal conformal transforma- 
tion z ^ z + cin.z"'^^ which is regular at the origin, 

0(0) ^ 0(0) + a_i[L_i0](O) + ao[Lo0](O), (34) 

with no further terms. 

^°It is possible to show that Eq is in general independent of a (in CFTs in which all 
the scaling dimensions are non- negative), so this is consistent. However, there is useful 
information in the L-dependence: if we normalise so that Eo{oo) = 0, then Eo(L) = 
— 7rc/24L, where c is the central charge of the CFT. 
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However, a more general conformal transformation of this kind is bound 
to have singularities elsewhere in the complex plane (including infinity). In 
order to understand the implication of this, let us consider a more general 
infinitesimal transformation in the form of a Laurent series, which may be 
singular at the origin: z z + Yl,^=-oo ^nZ"^^^ ■ No matter how small we 
choose the parameters a„ for n < — 1, for sufficiently small z this transfor- 
mation blows up. We deal with this by imagining cutting out a small half-disc 
of fixed radius around the origin, and, inside this, smoothing out the trans- 
formation. By definition it cannot be conformal in this region. The effect of 
this from the point of view of the correlation functions with other operators 
far from the origin is to generate new operators, called the descendents (p. 
Thus in general we have 

oo 

0(0) -.0(0)+ «nM](0) (35) 

n=— oo 

Notice that by dimensional analysis the descendent operators L„0 have scal- 
ing dimensions + (— n). 

We can keep doing this, by examining the properties of these descendent 
operators under further infinitesimal conformal transformations. We thus 
generate an enormous set of descendent operators Ln-^^Ln^ ■ ■ ■ (p, with scal- 
ing dimensions + This, in turn, implies that to each simple 
eigenstate of the strip transfer matrix H corresonds a whole tower of other 
eigenstates, with eigenvalues spaced by integer multiples of n/L.f^ 

^-'^These descendent operators do not necessarily satisfy the condition that [Ln(j>] = 
for n > 1, that is they are not primary and they do not satisfy (|2^) without adding some 
additional terms. In general, however, they are related to other descendent operators, with 
lower scaling dimension, are already constructed. For example, = 2[Lo4']- This 

algebraic structure is realised by thinking of the L„ as operators which act on the states 
of the transfer matrix. They are in fact the Fourier components Ln = e^'^™^^T{v)dv 
of a special operator called the stress-energy tensor. They satisfy the Virasoro algebra 
[Ln, Lm] — {n — m)Ln+m + ■^n{TT? — l)(5„^_m, which plays a very important role in the 
algebraic development of conformal field theory. However, we do not really need this 
formalism for the purposes of these lectures. 
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6 Null operators and differential equations. 



In general, with each scahng operator is associated an infinite number of 
descendents. However, they do not all have to be independent, and in many 
important applications, they are not. In that case, as we shall show, the 
correlation functions of cj) satisfy simple linear differential equations. 

Let us consider the case of interest, the bcc operator <f>f\a, with scaling 
dimension x{Q). Consider the partition function Zfa of an annulus <v < 
L, < u < W, with periodic boundary conditions in the li-direction, free 
boundary conditions f on v — and fixed boundary conditions a on v — L. 
As discussed above, we may write 

ZfaiQ) = TV eM-WHfa{L)) (36) 

j 

where q = e"'^^/^ and the sum is over all scaling operators, and their de- 
scendents, which may occur consistent with the given boundary conditions. 
The lowest is Xq — x{Q), corresponding to the operator itself. Let us 
examine the first few eigenstates above this, bearing in mind the fact that, 
as Q — > 1, Zfa{Q) — > 1. The first consequence is that, as already noted 
limQ^ix{Q) — 0. At Q — 1 this gives a contribution 1 to Z. All other 
contributions must therefore cancel among themselves in this limit. There is 
a (unique) operator L_i0j|a, with scaling dimension x{Q) + 1, which gives a 
contribution q to Z at Q — 1, which must be cancelled. The only way is by 
other primary operator (s) ip whose scaling dimensions x^{Q) — > 1 as Q — > 1. 
In fact, there are candidates for such operators: close to the bcc operator (f)f\a 
the Sq permutation symmetry of the Potts model is broken down to -S'q-i. 
The lattice operators 5s{r),b with h ^ a thus span a ((5 — 2) -dimensional space. 
We identify ■0 with the continuum limit of these operators, with degeneracy 
(Q-2). Thus 

Zf\a{.Q) = q""^^^ + q"^^^^^ + {Q- 2)g"'^(«) + • • • (37) 
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So far, so good.0 However, at the next level we have in principle the opera- 
tors L^2(pf\a, L'^i4>f\a (each with weight 1), and L^iip, with weight {Q — 2). 
Together these would give a contribution to Zaf. There are two possible 
resolutions: either there must be yet another new primary operator ip' with 
dimension x^'{Q) — > 2 as Q 1 and degeneracy — 1 + 0(Q — 1); or the 
two operators L_20/|a and L'^_i(t)f\a are not independent. In that case the 
combination L_20/|a — K,L'^i4'f\a (for some number k) corresponds to a null 
state, which does not contribute to the partition function. We shall assume 
that the latter solution is chosen. 

Let us now go back to the half-plane and examine the consequences of 
this for the four-point function 

G{Zi, Z2, Zs, Z4) = {(l)f\a{Zl)(l)a\f{z2)(l)f\a{z3)(j)a\f{z4)), (38) 

where the points Zj lie on the real axis. Consider the infinitesimal conformal 
transformation 



' = z + a{z) = z + a_2{z - ziy^ (39) 



z ^ z 

At the points Zj (2 < j < 4) this is regular, and the corresponding opera- 
tors transform simply according to (Pl). At zi however it is singular, and 
must be regularised in the manner described earlier. This then generates the 
new operator a^2[L-2(pf\a]{zi). We may now equate this to K[L'^i(f)f\a]{zi) = 
K{d/dziY(j)f\a]{zi), where k is an (as yet) unknown constant. Setting the 

Note that we are ignoring the possible contributions of operators with non-integer 
scahng dimensions at Q = 1. These are beheved not to arise, but in any case would not 
affect the argument. 

^•^This is the weakest point in the whole argument. However, there is further supporting 
evidence: at Q = 2 and at Q = 3 the same property is known to hold[|l^. In addition, 
the operator product expansion of two such operators is known, by the so-called fusion 
rules f^, t o generate another operator degenerate at level 3: this has been identified much 
earlier |l4| with the magnetisation operator at a free boundary, which, by duality, has the 
same scaling dimension as the bcc operator 4>a\b- 
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total variation with respect to a_2 to zero, we therefore find that G satisfies 
the differential equation 




x{Q) 



G{zi,Z2,Zs,z^) =0 (40) 



This is a linear partial differential equation. However, if we substitute 
in the form (p6D (which, recall, followed from the covariance under Mobius 
transformations), we find an ordinary second-order linear differential equa- 
tion for -F(?7) of a type very familiar to mathematical physicists. In fact, it 
is a hypergeometric equation, with regular singular points at (0,cxd,1). A 
lot of tedious algebra can be avoided by identifying the indices of this equa- 
tion, that is, the possible power law behaviour of the general solution near 
the singular points. For example, if F{ri) ~ as r/ — > 0, this implies that 
G ~ {zi2Z34)°'~'^^^'^\zi3Z2i)~°' in the limit where ^12^34 -C zi3Z24:- Substitut- 
ing into (pO]) and equating the coefficient of z'^2^^^^~'^ ^^ro gives 

K(a - 2x)(a - 2x - 1) + (a - 2x) - X = (41) 

Now we know that a particular solution is G, which in this limit should 
factorise: 

G ~ (0/|a(;^l)0a|/(^2))(0/|a(^3)0a|/(^4)) ~ ^r2''^'^^^34''^'^^ (42) 

This is because, in that limit, clusters which touch both segments (^1,^2) 
and (^3, Z4) are very rare. Therefore the indicial equation (^TJ) must have a 
solution with a = 0. This implies that 

3 

= 2{2x{Q) + 1) ^^^^ 

The other solution is then a = | + |a;(Q). 

A similar analysis may be made of the indicial equations at the other 
singular points. At 77 = 1 we also find (0, | + |a;((5)), and at 00, (— 4x(Q), | — 



25 



HQ))n 

The standard form of the hypergeometric equation is 



- r])F" + [c - (a + 6 + l)r/]F' - a6F = (44) 

for which the indices are (0,1 — c), (a, 6) and (0, c — a — 6) at (l,oo,0) 
respectively. This we have in our case 

a = -4x(Q), 6=|-|x(Q), c = § - |x(Q) (45) 

This equation has two independent solutions. The one which is regular at 
?7 = is proportional to 

^/ , N cib ri aia + l)h{h -\- 1) rf . „, 

.F.(a.t;c;,)^l + -l+ ^ |, + ■ ■ ■ (46) 

which converges for \r]\ < 1. There are in fact many different ways of writing 
down two independent solutions: a particularly useful pair are 

2Fi{a,b; a + b + 1 - c; 1 - rf) 

= 2i^i(-4x, |-|x;|-|x;l-r/) and (47) 

(1 - i]y''''\Fi{c - 6, c - a; c - a - 6 + 1; 1 - r/) 

= {1-11)1 + 1^2^1(1 - fx, 1 + fx; I + fx; 1 - r^) (48) 

Although we have been considering the four-point function of bcc operators 

corresponding to Zaa, the differential equation should hold equally well for 

that corresponding to Zab- In fact, the two independent solutions above are 

(apart from constants) just the solutions Faa and Fab corresponding to Zaa 

and Zab respectively. This is because rj 1 corresponds to Z2 z^. In 

the second case (ah), this introduces the non-trivial bcc operator (j)a\b at Z2, 

^"^The indices at 77 = oo are conventionally defined by F ~ 77"". There is a result that 
the 6 indices of such an equation must sum to 1. This explains the appearance of the 
fraction i when x{Q) = 0. 
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so cannot contain a term oc (1 — 77)° (incidentally, we see that this operator 
must have dimension | + ^x{Q)); while in the first case, the opposite is 
true. The coefficients may then found from the requirement that, as 77 — > 0, 
Faa ~ Fab 1, and the identity 

^/ , N T(c)T(c — a — b) . , 

2F1 a, h- c- 1 = ) ' ^ -( 49 

r(c — a)T[c — b) 

We therefore have the final results for the partition functions in the half- 
plane (normalised so that Zf = 1) 



Zaa = C^-(Q) ^;3 + 3^JU3 ^ fJ^ ,F,{-4x, 1 - |x; | - fo:; 1 - (50) 



,r(| + |x)r(i-|a;) 



r(i + x)r(| + |x) ,i + ^x 
" ^ r(| + |x)r(i + |x 



3 ' 3*^^ V3 ' 3- 

X 2Fi(| - fx, I + fx; I + fx; 1 - r^) (51) 

where C = {Z13Z24/ Zi2Z23Z3iZu). 

Actually, since we are interested in the limit Q ^ 1, that is x{Q) — > 0, 
these simplify greatly to the order needed:^ 

Zaa = ( '''''' y^^'^\l + %xhF,i-Ax,l;l;l-v)+Oix'){52) 
Zab = ^;^^|^(l-r/)l2Fi(i,|;|;l-r/)+0(x) (53) 

From d^), evaluated at Q = 1, now follows the main result for the crossing 
probability in the half-plane 



P((zi, Z2y, iz3, z,)) = 1 - ^^;j^|i, (l - r^) 32Fi(|, |; |; 1 - r]) (54) 
which may also be written 



P{{zi, Z2)\ (Z3, ^4)) = f(4jf(iy^3 2Fi(|, |; |; i]) (55) 



^^Notice that the ratio of gamma functions in Zaa simphfies to sin7r(| — \x)l sin7r(| + 
|a;), so is easy to expand in powers of x. 



27 



Although this follows from standard identities on hypergeometric functions, 
there is a more physical argument, based on duality in bond percolation: 
whenever there is a cluster of open bonds connecting {zi, Z2) to (2:3, Z4), there 
cannot be a cluster of open dual bonds connecting (^2, -23) to {z^, zi), and vice 
versa. Therefore there is either a cluster of open bonds connecting (zi, Z2) to 
(z3, 2:4), or a cluster of open dual bonds connecting (2:2, 23) to (2:4, zi). Thus 

P{zi, Z2] 2:3, Za) = 1- P{Z2, Z-i\ Zi, Zi), (56) 

or, in terms of the cross-ratio, 

P{r^) = l-P{l-r^). (57) 

The expression for the mean number of crossing clusters is slightly more 
difficult, because we need to differentiate with respect to Q, that is, with 
respect to x. First we note that Zaf and Zfa are equal to Zi2^ and z^l^ 
respectively. Thus the ratio of partition functions in (|1^) may written entirely 
in terms of 7] 

= (1 - r/)-2^(«)(l + i|a:(Q))2Fi(-4x(Q), i; |; 1 - rf) (58) 
to the order required. Thus 

- r(|)r(| + n)(i-r^)-' 
n^ir(|)r(i 

where the series comes from differentiating 2-^1 (a, • • •) with respect to its first 
argument a at a = 0. 

Of course this result depends on the (presently) unknown value of x'(l) = 
{dx/dQ)\Q=i. This may be founcQ by noting that, as — 0, it is ex- 
tremely unlikely that there is more than one crossing cluster, so that E[Ni)^ ~ 

In fact, the complete dependence of x{Q) has been conjectured using Coulomb gas 
methods. [r1 



E[N,] = x'{l) 



-2 Infl - 77) + ^ - 4 V ' 



(59) 
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1 

P ~ (r(|)/r(|)r(|))?73. Applying Stirling's formula to the coefficients of 

_4 

the series in (^) we see that they behave at large n like (r(|)/r(|))n 3, 
which implies that the singular part of the sum of the series behaves like 
(r(|)r(— i)/r(|))r73 as — > O. Equating this to the leading term in P then 
gives 

x'(l) = - (4r(|)r(-i))~' = -sin(47r/3)/(47r) = v^/(87r) (60) 

This result for x'(l) has a direct physical meaning. Consider once again 
the annular geometry, formed by sewing together the ends of a strip of width 
L and length W. The mean number of clusters crossing from one edge to 
the other is once again given by a formula like ([T^) , where now the partition 
functions may be represented in the trace form like (^). In the limit S> 1/ 
the result is particularly simple, since Zaa ~ Zjf — > 1, while Zaf = Zfa ~ 
exp{—7ix{Q)(W/L)). Thus in this limit, 

E[N,] ~ 2Trx'{l){W/L) = {V3/4:){W/L) (61) 

It is to be expected that, for large W, E[Nc] should be proportional to W, 
and therefore by scale invariance should go like (W/L), but the (irrational) 
universal coefficient is quite surprising. 



7 Other approaches. 

7.1 Crossing probability in the rectangle and modular 
invariance. 

An important special case is when the curve F is a rectangle, with the seg- 
ments 7i and 72 being opposite edges (see Fig. The conformal map- 
ping from the half-plane into the interior of the rectangle has the form of a 
Schwartz-Christoffel transformation 

""^"^^ Jo (l-t2)l/2(l_pi2)l/2 ^62) 
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which maps the points {—k ^,—1,1, A; ^)to the corners of the rectangle. The 
width and height of the rectangle are then 

/•I Hf 

W = 2 , = 2K{k^) (63) 

^0 ^(l-t2)(l_pt2) 
rl/k rjf 

L = . = K{1 - k') (64) 

•^1 ^{f^ -l){l-kH^) 



where K is the complete elliptic integral of the first kind. The aspect ratio 
r = W/L = 2K{k'^)/ K{1 — k'^). Clearly it is possible numerically to solve 
this for k for a given r, and compute 1] = ((1 — A;)/(l + /c))^ and thence 
the crossing probability P(r) via (|^). However, it turns out that there is 
a more elegant formula, in terms of modular forms, due to Kleban[|T3| (see 



also Ziff[|Ig].) Kleban and Zagier[T3] have recently derived this result, using 
duality and some basic properties of conformal field theory, with one rather 
simple assumption, and we now summarise their argument. 

From (H), P(r) = 1 — Zab, where Zab is the partition function for a strip of 
width L, with free boundary conditions on each edge, and length W. Think 
of evaluating this using the continuum version of the transfer matrix acting 
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horizontally, along the strip: 

^a(.= (a|e-^^//(^)|6) (65) 

where (a| and \b) are so-called boundary states, whose form we shall not need. 
Inserting a complete set of eigenstates into this matrix element: 

Zab = E{^\j)q''^^'{j\b) (66) 
j 

where q = e^'^^i'^/^) = e^'^^^ and the sum is over a complete set of scaling 
operators, consistent with the boundary conditions (//) on the edges, the 
lowest of which has x = 0. The sum may be organised into a sum over 
primary operators, and their descendents, so that P(r) has the form 

P(r) = 1 - E 1''^'QAq) (67) 

where tt// is the set of primary scaling dimensions consistent with the bound- 
ary conditions (//), and Qx{q) is a power series, starting with a term 0(g°), 
representing the contribution of all the descendents. Note that in this case no 
odd powers of q^^^ enter. This is a consequence of the reflection symmetry of 
the rectangle in the horizontal axis. Odd powers correspond to states which 
are odd under this operation, so do not couple to the states \a) and \b), which 
are even. On the other hand, we could also imagine computing Zab using the 
transfer matrix acting vertically. The result would have the form 

Zab = E r^'R^Q'^') (68) 

where q = e"^'^/'' and now odd powers of q^^'^ are expected to enter because 
Zab is not symmetric under reflections in the vertical axis. Note that the 
lowest value of x in TTab is the scaling dimension of 0a|b, which is strictly 
positive. 
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Finally, the function P{r) also satisfies (p7D, which means that P{r) = 
1 — P{l/r). Comparing with ( |68D we see that P(l/r) has the form of the 
series on the right hand side, and comparing this with ( p7| ) we see that in 
fact we may write 

P(r) = y: q'^'QM (69) 

XGTTab 

Now assume that there is only one term in this sum, corresponding to 
X = Xa\b- Let r = it, so that q = e^'^*'^. Then P is holomorphic in the 
upper-half r-plane. Define /(r) = dP/dr. The modular group is generated 
by the elements S : t ^ —l/r and T : r ^ r + 1, satisfying (ST)^ = 1. Now 
( |57D implies that under 5*, /(— 1/r) = —T'^f{T), while under T, /(r + 1) = 
^i-Kxji^^y Since {STY = 1, (e*'^'^)^ = —1, so that x = m/3 for some odd 
integer m. It follows that the function g = satisfies 1/r) = r^'^gij), 
g{T+l) = g{T), and that it has a series expansion in positive integers powers 
of q. It is called a cusp form of weight 12. It turns out that there is only 
one such function, up to an overall multiplicative constant, which is ri[qy^, 
where t] is the Dedekind function 
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7^ = gi/^^n(l-0 (70) 



n=l 

Thus 

27/3^2 



^ , ,r2 POO 

nr) = ^^^^/ V{^rrdr' (71) 
where the constant is fixed by the requirement that P{0) = 1, and various 



identities satisfied by ri{q). Remarkably, this formula is identical |T5| to the 
Schwartz-Christoffel transform of (^. 

7.2 Carleson's formula for an equilateral triangle, and 
Smirnov's proof. 

L. Carleson made the interesting observation that the crossing formula in an 
equilateral triangle has a much simpler form than that for a rectangle. This 
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Figure 7: Crossing cluster in the equilateral triangle. 



is because of a simple property of the hypergeometric differential equation 
arising from (^0]): because there is a solution with F = const., it follows 
that F'{ri) satisfies a simple first order equation which may be solved by 
quadrature. The result is that P{ri) may be written as 

But the occurrence of the exponent | means that this integral has just the 
form of the Schwartz-Christoffel conformal mapping from the upper half- 
plane to an equilateral triangle! 

Thus if we consider an equilateral triangle ABC (see Fig. |^), whose side 
has a unit length, and consider the crossing probability from anywhere on 
side AB to the segment XC of length x of the side BC, then Carleson's 
version of the crossing formula is simply 

P{AB,XC) = x (73) 
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Very recently, S. Smirnov]^ has provided a proof of this formula starting 
from site percolation on the triangular lattice. He considers the probability 
h{r) that a given point in the interior of the triangle is separated from the 
segment AC by at least one cluster spanning from AB to BC. Taking r = X 
on the boundary segment BC gives the required crossing probability. He 
shows that h{r) is, on the lattice, an approximate harmonic function, which, 
in the limit when the lattice spacing a — > 0, satisfies Laplace's equation with 
the boundary conditions that /i = on AC, h = 1 when r = B and the 
components of Vh at 60° to the boundary (i.e. parallel to AC) vanish along 
AB and BC. For the equilateral triangle, the solution of this is simple: h{r) 
is the distance of the point r from the edge AC, in units where the height 
of the triangle is 1. Of course, in general the solution of this boundary value 
problem is conformally invariant, so in the half-plane, one gets the result (ff^). 
Smirnov's arguments explain why the crossing formula (0) is the boundary 
value of an analytic function, but do not immediately give insight into the 
conformal field theory structure underlying this result. It seems to be special 
to Q = 1. 



7.3 The approach of Lawler, Schramm and Werner. 



Finally we discuss briefly another approach initiated by Schramm ||18||, which 
fits into a larger investigation into the intersection properties of Brownian 
walks by the above authors (see Ref. [0 and many references therein) The 
idea is to study the properties of a percolation hull, or cluster boundary, 
which grows from the real axis up into the upper half plane. Suppose that 
all the edges along the negative real axis are open, and those on the positive 
real axis are closed. There is an infinitely large cluster formed by the union of 
the links on the negative real axis and those clusters attached to it. Consider 
exploring the boundary of this cluster with a time-dependent process which 
makes one step along this boundary in unit time. (The process is almost 
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Figure 8: The path which explores the hull of a cluster attached to the 
negative real axis. 

symmetric under x — ^ —x because this path is also the boundary of a cluster 
of open dual links attached to the positive real axis.) The end of this path 
executes a non-crossing random walk in the upper half-plane, occasionally 
touching real axis (see Fig. H). This successively excludes a larger and larger 
region of the half-plane. Schramm argued that the continuum limit of this 
process may be defined in terms of the conformal mapping gt{z) which maps 
the complement of the excluded region onto the upper half-plane. He conjec- 
tured[^ that this function satisfies a kind of stochastic ordinary differential 



where go{z) = z and a{t) is a one-dimensional Brownian motion, i.e. a = ({t) 
where C is a gaussian noise satisfying ({t)({t') = K5{t — t'). This equation is 
defined up to the time that gt{z) hits a{t): we can think of this as the time 
Tz at which the point z enters the excluded region.^ 

^''The recent work of Smirnov|^ shows that percolation does indeed have a conformahy 
invariant continuum hmit. Uniqueness arguments then justify this conjecture a posteriori. 

^^It is useful to think about the case when a{t) = 0. The solution is then gt{z) = 
(z^ + 4:ty/^. The end of the path lies at 2i^/i at time t. The initial effect of the noise in 



equation 




(74) 
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The continuum problem is conformally invariant by construction. If we 
now consider the crossing probabihty between the intervals {—oo, —a) and 
(0, b), it may be related to the times T_a and Tb by 

P((-oo, -a), (0, b)) = Pr{T^a < n) (75) 

From the conjectured equation (0), with k = 6,0 it may be shown that the 
above is given by the crossing formula (0) with the cross-ratio rj = b/{a + b). 
Of course, many other related formulae also follow from this approach. 



8 Conclusions. 

The crossing formula @ was first conjectured in 1992. The original argument 

was based on (1) the mapping to the Q-state Potts model in the limit Q ^ 1 

and (2) ideas of conformal field theory which had been developed up to that 

time, which have been the main subject of these lectures. However, it is 

clear that percolation is somewhat special as in some sense the measure, 

suitably defined, is conformally invariant, rather than merely covariant as 

is believed to occur in most critical systems. This is linked to the fact that 

for percolation the central charge c vanishes.]^ This invariance has allowed 

mathematicians to formulate other more direct approaches to the problem. 

It remains to be seen whether these methods will be successful in computing 

other quantities, such as the mean number of crossing clusters (^, since, at 

least in our approach, this requires going away from Q = 1, and whether, 

in a more general context, they will shed light on the origin of conformal 

covariance in more general critical systems. 

a{t) is to give this an additional horizontal motion. 

^^Other values may be considered, which give generalised crossing formulas. In general, 
the process generated by ( |7^ ) is called Stochastic Loewner Evolution with parameter k 
(SLEk). 

^°In general the partition function of a finite system in a region size L with smooth 



metric and smooth boundaries behaves like Z ~ g^xi/e ■v\r];iere x is the Euler character.! 12 
Thus the measure is not strictly even scale invariant when c ^ 
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